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Measurement of the dark matter velocity anisotropy 
profile in galaxy clusters 
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Dark matter particles form halos that contribute the major part of the mass of galaxy clusters. The 
formation of these cosmological structures have been investigated both observationally and in nu- 
merical simulations, which have confirmed the existence of a universal mass profile [1]. However, 
the dynamic behaviour of dark matter in halos is not as well understood. We have used observa- 
tions of 16 equilibrated galaxy clusters to show that the random velocities of dark matter particles 
are larger on average along the radial direction than along the tangential, and that the magnitude 
of this velocity anisotropy is radially varying [2]. Our measurement implies that the collective 
behaviour of dark matter particles is fundamentally different from that of normal particles and the 
radial variation of the anisotropy velocity agrees with the predictions of numerical simulation [3]. 
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Dark matter is most often assumed to be a new type of particle which couples only very 
weakly to itself or other particle species [4]. Unlike a collisional gas, dark matter structures may 
therefore have a non-zero velocity anisotropy defined as jS = 1 — of /a, 2 , where a 2 and a 2 are the 
1 -dimensional tangential and radial velocity dispersions of the dark matter. Here we summarize the 
measurement of the radial velocity anisotropy profile in a sample of galaxy clusters, as presented 
in [2]. 

Most of the ordinary baryonic matter in galaxy clusters is found in the intracluster medium 
(ICM), a hot plasma emitting X-rays. We consider a sample of 16 relaxed galaxy clusters for 
which the radial ICM density and temperature profiles have been obtained from X-ray observations 
in earlier work. The clusters are selected to appear close to circular in projection, have smooth gas 
temperature and density profiles, and reconstructed total density profiles which are monotonically 
declining. The sample consists of eleven highly relaxed cool-core clusters at low redshift observed 
with XMM-Newton (A262, A496, A1795, A1837, A2052, A4059, Sersic 159-3, MKW3s, MKW9, 
NGC533, and 2A0335+096) [5, 6], and five intermediate redshift clusters observed with Chandra 
(RXJ1347.5, A1689, A2218, A1914, and A611) [7]. Data for A2052 and Sersic 159-3 were also 
used in an earlier analysis [8] where a constant velocity anisotropy was assumed. 

We determine the velocity anisotropy for the 16 clusters of our sample using the method out- 
lined in [2], and in almost all cases the anisotropy profile is radially increasing from close to zero 
in the center to about 0.5 or greater in the outer parts. Since the qualitative behaviour of the profiles 
are similar, we combine all our data into a single 'stacked' profile depending on radius in units of 
r 2500- Figure 1 shows the resulting velocity anisotropy profile which varies between 0.1 and 0.5 in 
the radial range 0.03-1. Or25oo. where ^2500 is the scale radius enclosing a mean density 2500 times 
greater than the critical density at the redshift of the cluster. Our measurement is in good agreement 
with the predictions of numerical simulations (see below) in the overlapping radial range. However, 
at large radii > r25oo> we overestimate the velocity anisotropy due to a lack of hydrostatic support, 
as demonstrated in [2], 

In order to calculate the velocity anisotropy, it is necessary to constrain the energy content of 
the dark matter. For this purpose we have assumed that the dark matter specific energy is equal to 
that of the ICM gas at the same radius up to a factor K, of order unity. This assumption is validated 
by the results of two independent numerical simulations of the formation of galaxy clusters in 
the ACDM cosmology, which both include dark matter and radiative gas physics. From each 
simulation we select the most relaxed clusters at redshift z = 0, which yields two samples of 67 
and 20 clusters that we refer to as the CLEF [9] and V06 [10-12] samples. For these two samples 
we calculate the two median K-profiles (Fig. 2) and we find that in both cases K is close to unity 
in the full radial range which is resolved, with about 10% scatter. Therefore, in our analysis we 
take K = 1.0 ±0.1. With this fc-profile, it is possible to calculate the radial velocity dispersion and 
the velocity anisotropy of the dark matter. This is done strictly using non-parametric Monte Carlo 
methods so that no modeling of the ICM or dark matter properties is involved. 

We have made a number of tests where we assume different fc-profiles, including both profiles 
that are greater or lesser than unity, and radially varying profiles which mimick the most extreme 
radial variations compatible with the CLEF and V06 fc-profiles. In particular, we check the effect 
of both AGN heating (k < 1) and radiative cooling (k* > 1) in the center of the clusters where 
simulations do not probe. In all cases the general behaviour is the same, i.e. a radially increasing 
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Figure 1: Radial profile of the stacked velocity anisotropy for our sample of 16 clusters. There are two 
slightly different ways of determining the velocity anisotropy. The error bars indicate the la percentiles of 
the underlying probability density. The hatched band shows the anisotropy profiles in the CLEF numerical 
simulation sample (see text), and the green line is the median anisotropy profile obtained when we apply our 
analysis to the CLEF sample. The correspondence between the reconstructed and actual velocity dispersions 
in the CLEF sample confirms the validity of our analysis. The leftmost vertical line indicates the innermost 
radius where the CLEF simulation resolves gravitational forces reliably, while the rightmost vertical line 
marks the onset of significant deviations from hydrostatic support in the CLEF sample. The horizontal 
dot-dashed line indicates the zero anisotropy of a gas. 




velocity anisotropy, although in some cases (k < 1) the velocity anisotropy can become zero at 
intermediate radii. See [2] for further details on these tests. 

Our result shows that the collective behaviour of dark matter particles is fundamentally dif- 
ferent from baryonic particles. Dark matter behaves as collisionless particles on the timescale of 
Gigayears, the dynamical timescale of galaxy clusters, and this implies an order-of-magnitude 
upper limit to the self-interaction cross-section per unit mass of o/ni < 1 cm 2 g , similar to what 
was found for the Bullet Cluster [13] and only slightly larger than the value proposed for self- 
interacting dark matter [14]. According to numerical simulations, the Galactic dark matter halo is 
also anisotropic, and this anisotropy will influence direct detection rates at the 5-10% level [15] 
and is potentially measurable in directional searches [16]. 
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Figure 2: Radial profile of K for the numerical simulation samples. The ratio of the specific energies is 
K = °dm/ ( J ? where /I = 0.6 is the mean molecular weight of the ICM gas particles. The vertical 
line indicates the largest radius of the observational data sample. 
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